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Abstract When two independently-prepared Bose-Einstein condensates are released 
from their corresponding traps, the absorbtion image of the overlapping clouds presents 
an interference pattern. Here we analyze a model introduced by Javanainen and Yoo 
[I], who considered two atomic condensates described by plane waves propagating in 
opposite directions. We present an analytical argument for the measurement-induced 
breaking of the relative phase symmetry in this system, demonstrating how the phase 
gets localized after a large enough number of detection events. 
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1 Introduction 

A long-standing problem in physics is to understand how relative phases are estab- 
lished between relatively independent systems that have never been in contact with 
each other. In Anderson's famous gedankenexperiment [2] two buckets of superfluid 
Helium are suddenly put in contact through a weak link: since there is no well-defined 
phase between the two condensate^], the question is then if a Josephson current will 
flow between them and why. More recently, experiments done at MIT [3] have estab- 
lished that clouds of Bose-Einstein condensates exhibit an interference pattern when 
released from the traps and imaged using a standard CCD camera. It soon became 
clear [TJH], that for standard single-shot absorbtion measurements the data cannot be 
explained simply by taking the average of the density operator; instead, what is in fact 
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1 It is important to notice that one should first cool the systems into the superfluid phase 
and only then connect them through the weak link. For otherwise, if the weak link is realized 
before, the system cools into the ground state, and a relative phase equal to zero will be 
established between the two components of the order parameter. 
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recorded are higher-order correlations. In a landmark paper by Javanainen and Joo 
P], a detection experiment on a system consiting of two counterpropagating beams 
is simulated numerically; it is found, surprisingly, that the result shows an interfer- 
ence pattern although the initial state is fragmented. In another important paper [4], a 
two modes/two detectors simplified model was introduced in analogy with the Hanbury 
Brown and Twiss experiment from optics. The essential feature captured by this model 
is that the detectors are such that they cannot distinguish between atoms prepared 
independently in two orthogonal modes. Other authors have extended the investiga- 
tion of the effects of measurement on various other types of many-body fragmented 
states involving spins [8] or rotating and attractive-interaction condensates [9]. One 
important prediction of these very elegant theories is that the interference fringes will 
be shifted randomly (corresponding to a random phase between the condensates) from 
one run to the other of the experiment. 

Although for the original MIT interference experiment [3] it can be argued that the 
randomness in the position of the fringes could also be caused by fluctuations in the 
trapping Hamiltonian [5] , a latter generation of experiments succeeded in getting better 
control over the external trapping parameters, as demonstrated by the reproducibility 
of the results when phase states are prepared and measured. Such interference experi- 
ments have been done in optical lattices [6] and with on-chip splitters [7] . The result is 
that indeed a random phase is detected for each realization of the experiment, leaving 
little doubt that the interpretation [TJ[4] is correct. Moreover, for some of these experi- 
ments the interaction energy can be neglected due to the low density of particles when 
the clouds corresponding to neighbouring sites overlap. 



2 Model 

In a real experiment, such as those mentioned above 3,6,7 , the atoms are released from 
a potential that confines them in the ground state, with an order parameter given by the 
Gross-Pitaevskii equation. In this initial state, the two clouds have mostly potential and 
interaction energy (as given for example by the Thomas- Fermi approximation); but, 
after the confining potential is switched off, this energy is rapidly transformed into 
kinetic energy. Depending on the particular setup, it can happen that at the time of 
the overlap of the two expanding clouds (i.e. of their time-dependent Gross-Pitaevskii 
wavefunctions) the interaction energy is completely negligible and therefore the atoms' 
motion is described by the standard Schrodinger equation. If this is not the case, the 
problem of taking the interaction into account is a challenging one [TO]. In this paper, 
we will neglect the interaction between atoms, and, to simplify the situation further 
[l] we will take the two atomic clouds as described by two plane waves of momenta 
±%k. A laser is then used to perform absorbtion (destructive) measurements, as in Fig. 
(fTJ). We note here that treating the problem without this simplification is possible by 
generalizing the results below, since approximate expressions for the order parameters 
of the expanding condensates are known. We work with periodic boundary conditions 
over a length L; therefore k = (2n/L)l, where I is an integer number. 
The initial state of the system is a fragmented or Fock state, 

| (N/2) k , (N/2)_ k ) = jJj^ b{ N/2 P_T | 0). (1) 
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Fig. 1 Schematic of the experiment: atoms are released from left and right with wavevectors 
k and —k repectively and a source S {e.g. a laser) is used to take images of the overlapping 
region by sending photons to a detector D {e.g. a CCD camera). 



The Hamiltonian of the system has only the term corresponding to kinetic energy, 

i dx ^(^ 2 Mx) = ¥-Y. ^ ( 2 ) 
J q 

so the atoms are not aware of each other's existence and as a result the state of the 
system at a later time t, satisfying the Schrodinger equation ih(d/dt) | W) = H | 
will be 

9) = e - lhk2m / 2m | (N/2) k , (N/2)_ k ), (3) 

where the phase factor is irrelevant. At the moment of detection, the state of the system 
is therefore a fragmented state. 

The field operator can then be expanded in the corresponding mode operators, 

p iqx 

and, taking into account that only k and — k are occupied, we can truncate the expan- 
sion to b k and b_ fe only, {j)(x) » (1 / '\fL)(b k e lkx + b_ k e~ lkx ). 

Fragmented states do not seem at first sight too exciting candidates for interference: 
consider | 9) = | {N/2) k , (N/2)_ k ) in the limit of a large number of particles N 1. 
Then the probability of finding a particle at a point x is the same for every point so 
the average density does not present any spatial structure: 

| 4> + (x)4<(x) | 9) = N/L. (5) 
Contrast that with the case of a phase state, 

| <p) N = _^L=(6+ e ^/ 2 + bt k e-^ /2 ) N | 0), (6) 



where the density 



N {cp | iP+{x)iP{x) | ip) N = — (1 + cos(2fcc + ip)). (7) 
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oscillates in space with a period given by the relative momentum 2Tik between atoms 
with wavevectors ±fc. However, for a fragmented state something interesting seem to 
happen with the second order (density-density) correlation function, 

;y2 -I 

{9 | i> + (x 1 )i, + (x 2 )^{x 2 )i,{x 1 ) \W) = —(! + - cos2fe(x 1 - x 2 )), (8) 

which reveals an interference structure. 

The same results can be obtained if the calculation is performed in the phase repre- 
sentation: a fragmented state can be expanded in phase states | <p) N = (2 N N\)' 1/2 [b+ e iip/2 + 
b+ k e-' lip / 2 ) N | 0), in the form g] 

Now, using the orthogonality relation jy{(p \ if')]\f = 2y^2ir/N5(ip — ip'), we can easily 
recover the same expressions for the correlation functions Eqs. (|5I8|I for the fragmented 
states with the additional benefit of the insight that integration over phases washes out 
the interference fringes in the density Eq. (0 for phase states. What does this mean? 
The interference pattern in Eq. ([8} is obtained by removing an atom at one position 
and adding another one at a different position (correlations functions of this type do 
indeed test the response of the system at one point given that an operator is applied at 
another point). This suggests that, in order to predict the result of such an experiment, 
we have to carefully consider the effect of the measurement process itself one atom at 
a time, which we do in the next section. 

To make the connection with other insights into this problem, we first look at the 
combination of operators in Eq. ^ which is responsible for the interference term, 
and we find the operator A = btb_f. introduced in [11], which can be regarded as 
the quantum observable giving the amplitude of the fringes; the square amplitude 
of the interference term in Eq. (J5J is obtained from averaging the operator A + A = 
—kk^ 1 " name ly I ^ifc^fc ^fc^-fe I ^) ~ (N/2) 2 . Second, we notice that the two- 
detector model presented in [4] can be understood as a particular case of the model 
above, namely that in which detection can be realized only at positions x+ = and 
X- — L/41 (kx~ = tt/2), leading to the "detector" operators \/L$(x+) = b^. + b_f. 
and —iyL'if)(x-) — b^ — &_fc. 

3 Measurement 

What happens during the measurement of such a state? We divide the length L into 
M q > 1 intervals of length L/M each, and we consider n atoms detected, n\ atoms 
in the first interval n 2 atoms in the second and in the last interval. Following the 
quantum- optics theory of detection, the many-body state after n = n\ + n 2 + .... + tim 
atoms are detected becomes 

M 

n^) ni i#>jv= 

i=l 

(i) V^rhr) L 2?n cosn "^+^ 2 )i^-«- (10) 



5 



The statement that we want to prove is that sequences of rii detections at positions 
Xi occur with maximum probability if the Born probability rule m = 2{n/M) cos 2 {kx,i + 
<p/2), is satisfied, where tp is an arbitrary angle. We proceed by mathematical induction: 
given that this statement is true for n detections, we want to prove that the next 
detection event will tend to occur with maximum probability such that the law above 
is satisfied. As the number of detection events increases, the next detection events tend 
to "reinforce" this law, in the sense that the probability for events that do not satisfy 
the Born probability distribution becomes exponentially small. The tendency of phase 
localization increases exponentially from one detection to the other. 

We start with the first two detection events and we go back to understanding the 
structure of Eq. ([8]). The probability of detecting a particle at position Xi with a system 
in a fragmented state is proportional to the first-order correlation which is flat. That is, 
all points Xi are equally probable, P[xi) = 1/M. Not so for the next detection event, 
for which the probability will be proportional to the second order correlation function, 

P ( x 3 I x i) = Jj (l + \ cos{2kxj - 2kx t )^ . (11) 

Thus, a single detection at Xj sets up a phase 2kxi which will make the next detection 
event more probable at points Xj where cos(2kxj — - 2kx^) = (1/2)(1 + cos 2 (kxj — 
kxi)) is large. This also means that after the first detection event the relative phase 
of the system is no longer undefined, which can be checked directly for example by 
calculating the fluctuations of the phase operator [12] in the state ip(xi) \ N/2, N/2) . 
A more elegant way is to employ Bayes' theorem [131114] and attempt to calculate the 
conditional probability that the relative phase takes a certain value tp given a detection 
event at x.^, P(tp \ Xi): the probability of obtaining a count at Xi is P(xi) — 1/M, 
and the probability of measuring a certain phase tp on the fragmented state is also 
equally distributed, P{tp) — l/2rr. If the system would have been prepared in a phase 
state | tp) N , the probability to count a particle at position xi had been P(xi \ tp) — 
M _1 [l + cos(2A;Xj + tp)]. Bayes' theorem allow us to change our past description of a 
system given that further information is available: in our case, it relates the probabilities 
above as 

P(tp | x % ) = JyprP{xi | tp) = + cos(2fcx. t + tp)], (12) 
P(Xi) 2-k 

showing indeed that the most probable phases are those that get close to 2kx.i. 

We now look at what happens if the statement above is valid up to the n'th detec- 
tion event. This means that we identified an angle tp n so that 2(n/M) cos 2 (kxi + tp n /2) 
closely resembles the histogram n^a^. 

We first notice that this results in an extremum for the probability amplitude: 
M 

= TT cos™' ( kx i + ¥>/ 2 ) = 
dtp J- J- 

i=l 

M 

- rii sin(fcj + tp/2) coa ni ' 1 {kx i + tp/2) TT cos™^ (kxj + tp/2). (13) 

It is immediate to see that this equation is satisfied for rii = 2{n/M) cos 2 {kxi + tp n /2); 
indeed plugging this expression into Eq. (]13[) results in 

d M 

— T\cos ni (kxi + <p/2) | p =£ n = 
aV i=i 
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n 
~M 



/ M \ M 

sin(2fc i a; l + 

/ i=l 



! J cos ni (kxi + ipn/2) = 0, 



(14) 



since the first sum vanishes. We note that our solution gives the correct normalization 



condition X]i=i n i ~ n ' an< ^ that nj/n = (L/M)(y/2/Lcos(kxi + ip n /2)) 2 is precisely 
Born's probability law corresponding to a wavefunction \J~2j~L cos(kxi + <p n /2). After 
n detections the amplitude becomes peaked around a certain value Cp n of the relative 
phase, which satisfies as closely n, = 2(n/M) cos 2 (kxi + <p n /2); further absorbtion 
events will make this amplitude even more peaked. The final result is that after a large 
enough number of detection events, the initial fragmented wavefunction collapses onto 
one of the phase states. (This result is rather general; it does not depend on the specific 
form of the wavefunction, as long as it is normalized [ID].) To convince ourselves that 
this is indeed the case, consider the square of the amplitude probability 



M 



(15) 



Since by definition this quantity is positive, we can take the In and expand it in a 
Taylor series. As expected, the first order derrivative is zero at <p = <p, indicating an 
extremum, 



dtp 



m /(^) \tp=<p n = 



M 
i=l 



= 0, 



while in the second order we get 



d 



— In f(<p) \ v ^ n = -n. 
This means that we have a indeed a maximum at tp and 



f(<p) = /(v)e-' 



(16) 



(17) 



(18) 



The essential thing to notice is that the variance of this normal distribution function 
is simply n^ 1 , increasing with the number of detection events but independent on 
any particular outcome (i.e. how many detection events n% we had at position a;,-). 
After only a few detection events this will result in a sharply peaked distribution 
around the value ip (which depends on the distribution {n^Xj}). Indeed, it has been 
demonstrated experimentally that extracting a small fraction of the total number of 
atoms is enough to establish a relative phase [15] • Assuming TV — n >> 1, and using 
the quasi-orthonormalization relation for the phase states, we obtain the normalized 
wavefunction in the form 



' 4 - Jo 



dipe 



■±n((p-ip n ) 2 



v) 



N- 



(19) 



We now take the n-times collapsed wavefunction | gA")) resulting after a sequence 
of {xj,rij} detections and ask what is the probability of detecting a particle at a 
position x. This quantity will be proportional with (fl^ n ) | (xi)'ip(xi) | W^- 71 '). Using 
the properties of Gaussian integrals, one can easily calculate 



4nN{N - n) 



L 



1 + e 



~l/2n CQS (2kxi + (p) 



(20) 
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resulting in a conditional probability 

(21) 



P(xi | {xj,nj}) = 1 + e 1 ^ 2 ™ cos(2fexi + <p n ) 



(or, equivalently, a conditional probability density L 1 1 + e 1 / 2 ™ cos(2fcri + (p) ). 
From this expression it follows that events which fall far from the Born probability 
distribution [-v/2/Z cos(fcrj + tp n /2)] 2 are very unlikely to occur. Consider for example 
a detection at a minimum of this probability, say Xi — (n — (p n )/2k; the probability of 
this event is 

P((n - p n )/2k | { Xj , nj }) = ^(1 - e- 1/2 ") "3 X 0. (22) 

Another way to look at this effect is to think in terms of probabilities of detecting a 
phase tp: from Eq. (|19|l above, it is obvious that the probability of detecting a phase tp 
is a Gaussian centered at <p n and of variance 



P(<p) = il*l e - n te-9~y. (23) 

We now try to get an estimate for the behavior of the optimal phase tp as the number of 
detection events increases. To do we that we evaluate the difference between the optimal 
phases for two consecutive detection events, <p n +i — <Pn- Suppose that a detection event 
occurs at xf, then we find the correction made to the optimal value tp by applying 
the operator %[> + (x^ to the state | tp^) and, as before, finding the extremum of the 
quantity cos(kxi + tp/2) exp[—n(tp—tp n )/2] appearing under integration over tp. We have 
also checked numerically that for n> 1 this quantity has a unique, sharp maximum 
value in the interval [0, 2ir], which is a solution of the equation 



Vn+l — <Pn 



.Itan(fcx i + ^tl). (24) 



An analysis of this expression shows that the new value (p n +l is within 1/n of the old 
value ip n for events that have a finite probability (i.e. for which [y // 2/L cos(kxi + tp/2)] 2 
is not too close to zero). This result is intuitively clear: if a detection event occurs at 
a maximum of [W2/L cos(fcxj + tp/2)] 2 , the phase is not modified much. More change 
could happen if some improbable event ([-^/2/L cos(fcrj + p/2)] 2 small) occurs, and 
this could produce significant changes in the phase only if n is not too large (in other 
words, the value of the phase is decided in the beginning: the more detection events 
accumulate, the less is the influence of an improbable event, even if it occurs). As a 
result, the evolution of the phase tp n is given by a discrete stochastic process, with 
steps typically of the order of 1/n given by Eq. (I24p and probabilities Eq. (|2ip . leading 
very fast to localization. 

In the case of the particular model of two-modes/two-detectors [4] we obtain for the 
conditional probabilities Eq. (|21|l that an atom is counted in the "+" or "-" detector 
given a previous sequence of n+, n_ detections (n = n+ + n_) as 

P(±|r I+ ,n_) = i(l±e- 1 / 2 ™), (25) 

and the equivalent of the equation Eq. (I24[) is 



Pn+l ~<Pn = T~ [tan((p„ + i/2)] ±1 . (26) 
n 
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The same phenomenon occurs here as well for relatively large n> 1: the next value 
fn+l will fall within 1/n of the previous value, approximating better and better the 
final value of the phase. To give a numerical example, let us suppose that at n = 500 
we have obtained ^500 = 0. A detection in the + channel will only "confirm" the value 
zero for the relative phase, since the equation <^soi = —2 x 10~ 3 x tan(<^5oi/2) has 
V?50l = as solution. This value can be changed though by a detection in the other 
channel: indeed, solving numerically the equation 1 = 500 x ^501 x tan(<^5oi/2) results 
in <f5oi = 0.063, a small correction but worth considering. However, the probability 
for a detection in the - channel is only P{—) = 5 x 10~ 4 . Clearly, to get a significant 
change in the value of the relative phase one would need a series of such events, which 
is very improbable. 



4 Conclusions 



We have presented an argument based on mathematical induction for the localization of 
phase under absorbtion measurement of a Fock state consisting of counterpropagating 
atomic plane waves. We show that, as the number of atoms detected increases, an 
arbitary phase builds up. 
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